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By numerical simulations, we investigate the onset of synchronization of networked phase oscillators under 
two different weighting schemes. In scheme-I, the link weights are correlated to the product of the degrees 
of the connected nodes, so this kind of networks is named as the weight-degree correlated (WDC) network. 
In scheme-II, the link weights are randomly assigned to each link regardless of the node degrees, so this kind 
of networks is named as the weight-degree uncorrelated (WDU) network. Interestingly, it is found that by 
increasing a parameter that governs the weight distribution, the onset of synchronization in WDC network is 
monotonically enhanced, while in WDU network there is a reverse in the synchronization performance. We 
investigate this phenomenon from the viewpoint of gradient network, and explain the contrary roles of coupling 
gradient on network synchronization: gradient promotes synchronization in WDC network, while deteriorates 
synchronization in WDU network. The findings highlight the fact that, besides the link weight, the correlation 
between the weight and node degree is also important to the network dynamics. 

PACS numbers: 05.45.Xt, 89.75.Hc 



In the past decades, complex networks have became an 
established framework to understand the behavior of a 
large variety of complex systems, ranging from biology 
to social sciences Ql!-13]- In network studies, one impor- 
tant issue is to investigate the influence of the network 
topology on the dynamics, e.g., the synchronization behav- 
iors of coupled oscillators |4]. For many realistic systems, 
besides the degree, the weight of the network links also 
presents the heterogeneous distribution, i.e., the weighted 
networks. Depending on the correlation between the node 
degree and the link weight, the weighted systems can 
be roughly divided into two types: weight-degree corre- 
lated (WDC) and weight-degree uncorrelated (WDU) net- 
works [3]. While the collective dynamics taking place on 
weighted networks have been extensively studied Hfl, but 
so far, to the best of our knowledge, the impact of weight- 
degree correlation on the collective dynamics has not been 
addressed in literature. In this paper, we consider the ef- 
fect of weight-degree correlation on the onset of synchro- 
nization of a generalized Kuramoto model with fixed total 
coupling cost. Interestingly, it is found that, by increasing 
the weight parameter that characterizes the distribution 
of links weights, the synchronization is monotonically en- 
hanced in WDC network, while the synchronization could 
be deteriorated firstly and then enhanced in WDU net- 
work. Moreover, we explain qualitatively the fundamental 
mechanism from the viewpoint of gradient network. Our 
findings may be helpful to the design and optimization of 
realistic dynamical networks. 
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I. INTRODUCTION 

Due to its broad applications, synchronization of coupled 
phase oscillators has been extensively studied in the past 
decades |Hl 0] . Recently, with the blooming of network sci- 
ence, this study has been extended to the situation of complex 
network structures, i.e., the generalized Kuramoto model O- 
Il5n . Different to the traditional studies of regular networks, 
in the generalized Kuramoto model more attention has been 
paid to the interplay between the network topology and dy- 
namics. Considering the fact that links in realistic networks 
are generally of different strength lU^-HH, more recently the 
study of onset of synchronization in weighted complex net- 
works has arisen certain interest, where a number of new phe- 
nomena have been identified |20->2^1. 

For weighted networks in reality, a typical feature is that 
the weights of the network links possess a large variation. In 
particular, it is found that in many realistic networks the link 
weights roughly follow power-law distributions, say, for in- 
stance, the airport network ll30ll . the scientific collaboration 
network ll3lll . and the mobile communication network ll32ll . 
to name just a few. In the previous studies of weighted net- 
works, the heterogeneity of link weights has been attributed 
to the heterogenous network structure, e.g., the heterogenous 
degree distribution, and the weight of a link in general is set 
to be proportional to the product of the degrees of the nodes it 
connects. For this feature, we call this type of network as the 
weight-degree correlated (WDC) network, which is often ob- 
served in the biological and technological systems, e. g., th e E. 
Coli metabolic network l33ll and the airport network fl30l l33l . 
Lately, with the progress of network science, it is also found 
that some realistic networks, while the link weights also pos- 
sess a power-law distribution, the link weight and node degree 
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are not closely correlated, e.g., the scientific collaboration and 
email networks l30l[34ll . More specifically, in this type of net- 
works the link weights seem to be randomly chosen within a 
certain range, regardless of the node degrees. For this prop- 
erty, we call this type of network as the weight-degree uncor- 
rected (WDU) network. Our main task in the present paper is 
to investigate the dynamical behaviors of these two different 
types of weighted networks, with the purpose of exploring the 
effect of weight-degree correlation on the network dynamics. 

To study the effect of weight-degree correlation on network 
dynamics, we will employ the model of networked phase os- 
cillators, i.e., the generalized Kuramoto model, and study the 
onset of synchronization between the following two weight- 
ing schemes. In scheme-I, the coupling weight between a pair 
of directly connected nodes in the network is set to be pro- 
portional to the degrees of these two nodes, so as to generate 
a WDC-type network. In scheme-II, the weights of the net- 
work links also have a power-law distribution, but the weights 
are randomly distributed among the links, regardless of the 
node degrees. This gives rise to a WDU-type network. Our 
main finding is that, by tuning a weight parameter (to be de- 
fined later) that governs the heterogeneity of the link weights, 
the two types of networks have very different synchronization 
responses. In particular, in WDC-type network, with the in- 
crease of the weight parameter, the network synchronization 
is found to be monotonically enhanced; while in WDU-type 
network, it is found that the increase of the weight parameter 
may either suppress or enhance the network synchronization. 
We analyze the different synchronization responses of the two 
networks from the viewpoint of gradient network Il25[l35ll . and 
get the general conclusion that synchronization is enhanced 
by coupling gradient in WDC network, while is suppressed in 
WDU network. These findings highlight that, besides the link 
weight, the correlation between the link weight and the node 
degree should be also taken into account when analyzing the 
network dynamics. 

The rest of this paper is organized as follows. In Sec. |TTJ 
we will introduce our model of coupled phase oscillators, as 
well as the two types of weighting schemes. In Sec. [TTH we 
will show the transition paths from non-synchronous to syn- 
chronous states for the two types of networks, in which the 
critical coupling strength characterizing onset of synchroniza- 
tion will be defined. In Sec. [IV] we will study in detail the ef- 
fect of weight-degree correlation on the network synchroniza- 
tion, in which the reversed synchronization in WDU network 
will be highlighted. Also, based on the concept of gradient 
complex network, we will give an analysis to the synchro- 
nization performance of the two types of networks. In Sec. [V] 
we will give our discussion and conclusion. 



H. THE MODEL 

In our study, we employ the following two typical mod- 
els for the network structure: the small-world and scale-free 
networks. The small-world network is constructed by the 
method introduced in Ref. 03611 . in which each node has 2m 
links and each link is rewired randomly with a probability 



p r . The scale-free network is generated by the standard BA 
growth model ll37ll . which has average degree (k) = 2m, 
and the node degrees roughly follow a power-law distribution, 
P{k) ~ fc-T, with 7 = 3. 

For illustration, we set the weights in the two types of net- 
works as follows. For the WDC network, motivated by the 
empirical observation in metabolic networks [33] and airport 
networks BOl 13311 . we set the weight, c mn , between the pair 
of directly connected nodes, m and n, to be 

where k m and k n are the degrees of the nodes m and n, re- 
spectively, and a G K is a tunable parameter characteriz- 
ing the weight distribution, i.e. the weight parameter. For 
the WDU network, motivated by the empirical observation in 
scientific collaboration networks ll30ll and email net works ll34ll . 
the weight of each link is randomly chosen within the range 
[cmin, Cmax], while keeping the weight distribution to be hav- 
ing the power-law scaling, 

P(c) ~ c-P. (2) 

In WDU network, the weight parameter is f3, which, according 
to the empirical observations, is defined to be positive. 

Arranging each node in the network with a phase oscillator 
and regarding the weight as the coupling strength between the 
connected nodes, we now can study the dynamical behaviors 
of the two types of networks constructed above. Our model of 
coupled phase oscillators reads 

£ N 
n— 1 

with N the network size, and m,n — 1, . . . , N the node 
indices. <fi m and oj m represent, respectively, the phase state 
and the intrinsic frequency of the mth oscillator (node). The 
weights of the network links are represented by the matrix 
{c TO n}> with c m „ = if m and n are not directly con- 
nected, e is a uniform coupling strength which, for the sake 
of convenience, is normalized by the node weight intensity 
s m = Sn=i c mn- Our main task is to investigate how the 
weighting schemes, i.e. weight-degree correlated or non- 
correlated, will affect the network dynamics. 

In next simulations, we fix the network size as N = 6400 
and the average degree as (fc) = 2m = 6. In constructing 
the small-world network, we set the rewiring probability p r as 
0.2. For each network example, the intrinsic frequency, {tOi}, 
of the network oscillators is randomly chosen from the Gaus- 
sian distribution g(uj) = (2-Ka 2 )^ 1 ^ 2 exp(— o; 2 /2ct 2 ), which 
has zero mean value and the unit variance (a 2 = 1). Also, the 
initial states of the oscillators are randomly chosen from the 
interval [0, 2tt). With these settings, we then set a value for 
the uniform coupling strength, e, and evolve the network ac- 
cording to Eq. (0. The equations are integrated by the fourth- 
order Runge-Kutta method and the time step is At = 0.01. 
After discarding a transient process of 2000 time steps, we 
start to measure the collective behavior of the network, which 
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FIG. 1 : (Color online) Under the WDC scheme, the variation of the 
synchronization order parameter, R 2 , as a function of the coupling 
strength, e, for different weight parameters, a. The network is gen- 
erated by the standard BA model (scale-free network), which has size 
N = 6400 and average degree (k) — 6. Onset of synchronization is 
identified as the coupling strength where R 2 starts to increase from 
2 x 10" 2 . The error bars are estimated by the standard deviation. 
Each data is averaged over 20 network runs. 
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is characterized by the network order parameter [9 , 1 3 J 



R 



Here (. . .) represents that the result is averaged over a pe- 
riod of 2000 time steps, and [. . .] represents the average of 
20 network runs. d n is the total incoming coupling strength 
of the node n, which, due to the normalization operation in 
the model, is unit. By scanning the coupling strength e, we 
can monitor the behavior of R, so as to find the critical cou- 
pling strength, e c , where the onset of network synchronization 
occurs. 



III. ONSET OF SYNCHRONIZATION IN WEIGHTED 
COMPLEX NETWORKS 

We start by demonstrating the typical path from non- 
synchronous to synchronous states in non-weighted complex 
networks, i.e. the case of a = in Eq. (fTJ, so as to define 
the onset of network synchronization. As plot in Fig. Q](the 
curve for a = 0), when the coupling strength is small, e.g. 
e « 0, the trajectories of the oscillators are almost uncorre- 
lated, and we roughly have R w 0{\/\N) for the order pa- 
rameter. This value of R keeps almost unchanged till a critical 
coupling strength, e c , is met. After that, as e increases from 
e c , the value of R will be gradually increased and, when e is 
large enough, it reaches the limit R — 1 (the state of complete 
network synchronization). The value e c is then defined as the 
onset point for synchronization, which, for the case of non- 
weighted scale-free network, is about 2 (Fig. [T). In numerical 
simulations, we define the onset of network synchronization 
as the point where R 2 > 2 x 10~ 2 . In what follows, we will 
study how the value of e c is changed by the weight parame- 
ters, a and (3, in the two types of networks. 




(c) 



FIG. 2: (Color online) Under the WDU scheme, the variations of 
the synchronization order parameter, R 2 , as a function of the cou- 
pling strength, e, for scale-free networks (a) (c) and small-world net- 
works (b). As a reference, the synchronization of non-weighted net- 
work (Binary) is also plotted in each subplot. All the networks have 
size = 6400 and average degree 6. In constructing the WDU 
network, each link is arranged a weight chosen randomly from the 
range [1, 100] for (a) and (b), and from the range [1, 100000] for (c). 
The rewiring probability in generating the small-world networks is 
p r = 0.2. The error bars are estimated by the standard deviation, 
and each data is averaged over 20 network runs. 



We first check the influence of the weight parameter, a, on 
WDC-type networks. The results are carried for scale-free 
networks and are plotted in Fig. Q] It is shown that, as the 
weight parameter increases from —3 to 5, the network syn- 
chronization is monotonically enhanced. In particular, the 
point of onset of synchronization, i.e. the vale of e c , is found 
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to be gradually shifted to the small values. For instance, when 
a = — 3, the critical coupling strength is about 7, while this 
value is decreased to about 0.9 when a — 3. Another observa- 
tion in Fig. Q]is that, in the regime of a < —1, a small change 
of the weight parameter could induce a significant change to 
the network synchronization; while in the regime of a > —1 
the similar change in a has relatively small influence. For 
instance, by increasing a from —3 to —2, the critical cou- 
pling strength is decreased from 7 to 4.1; while it is almost 
unchanged when a is increased from 3 to 5 (see Fig. [TJ. 

We next investigate the influence of the weight parameter 
on WDU-type networks. In Fig. |2ja) and (c), we plot the 
variation of the order parameter as a function of the coupling 
strength for the scale-free network. It is found that, similar 
to the case of WDC networks, the network synchronization 
is also influenced by the weight parameter. The new finding 
is that, in WDU-type networks the influence of the weight pa- 
rameter on network synchronization is non-monotonic. In par- 
ticular, it can been seen from Fig. [2 a) and (c), as f3 increases 
from 0.8 to 1.4 the network synchronization is deteriorated; 
while increasing /3 from 1.4 to 2.4, the network synchroniza- 
tion is enhanced. Speaking alternatively, as the weight param- 
eter varies, there is a crossover in the synchronization perfor- 
mance for WDU network. Under the same weighting scheme 
(WDU network), the crossover of network synchronization is 
also observed in the model of small-world networks, which 
can be seen from Fig. EJb). In particular, the max weight c max 
has no effect on the onset of synchronization on scale-free net- 
works and small world networks (only give the examples on 
scale-free networks), qualitatively. 

The results of Figs. Q] and [2] provide us with the follow- 
ing information: 1) the weight parameters have important ef- 
fects on the network synchronization, and 2) the effects of the 
weight parameters on network synchronization are dependent 
on the weight-degree correlation. To explore these findings in 
depth, we go on to study the synchronization mechanisms of 
the networks. 



hanced. The influence of gradient on the onset of synchroniza- 
tion of coupled phase oscillators has been also investigated 
lfl3tl . where an important finding is that, besides the gradient 
strength, the gradient direction also plays a key role in influ- 
encing the network synchronization. In particular, in Ref. lfl3ll 
it is demonstrated that when the gradient direction is pointing 
from the larger-degree to smaller-degree nodes on each link, 
the increase of the gradient strength will enhance the synchro- 
nization monotonically; while if it is in the opposite direction, 
the increase of the gradient strength will suppress the synchro- 
nization. Noticing of these effects of the weight gradient and 
gradient direction on the network synchronization, we in the 
following explore the properties of the gradients in our mod- 
els of WDC and WDU networks, with a hope to explain the 
numerical observations in Figs. [T]and|2] 

Following the method of Ref. [25], we construct the gradi- 
ent networks in our models, as follows. The Eq. (|3) can be 
rewritten as 



N 

£ ^ 9mn sin( 

n=l 



(5) 



with g mn = c mn /s m the normalized coupling strength, which 
captures the real coupling strength that node m receives from 
node n. The reason for the normalized coupling strength is 
simply for convenience, since in this way the total network 
coupling strength will be keeping as constant and is not chang- 
ing with the weight parameters. That is to say, in our models 
the change of network synchronization is caused solely by the 
weight distribution, instead of the change of the coupling cost. 

The coupling gradient, Ag mn , flowing from node n to node 
to under the condition k n > k m is written as 
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IV. THE INFLUENCE OF WEIGHT-DEGREE 
CORRELATION ON NETWORK SYNCHRONIZATION 

In studying the dynamics of directed and weighted complex 
networks, an effective method is to analyze the properties of 
the weight gradient ll25l l26ll . For a directed network, the pair 
of weights on a link in general is different from each other, 
which causes a gradient (bias) among the connected nodes. 
The idea behind gradient analysis is to extract the asymmet- 
ric part (gradient) on each link, and analyze its functions to 
the network dynamics separately. In this way, the original 
network can be virtually regarded as a superposition of two 
subnetworks: one constituted by the symmetric links and the 
other one constituted by gradient links. In Ref. ll25ll . the later 
is also called the gradient complex network. 

Gradient network has been employed in analyzing the net- 
work synchronization in literature. In Ref. H25II it is shown 
that, with the increase of the gradient strength, the complete 
synchronization in scale-free network can be significantly en- 



where c mn = c nm is the weight before the normalization op- 
eration, and s is the node weight intensity defined in Sec. 
II. Extracting only the gradient element on each of the net- 
work links, together with the network nodes, we can construct 
the gradient network. In our study, the positive direction of 
the gradient is defined as pointing from the larger-degree to 
smaller-degree nodes on each link. With this definition, we 
can measure the gradient at the network level, i.e. the aver- 
aged network gradient, which is 
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where (. . .) means that the value is averaged over the network 
links (totally there are 2M — N(k) links in the network). 
Please note that the summation in Eq. (7) is under the condi- 
tion of k m < k n , which is based on the direction of coupling 
gradient defined above. Our next job is to analyze how the 
tuning of the weight parameter will change the gradient net- 
work and, in turn, affect the network synchronization. 



5 



A. WDC-type Networks 

We first analyze the synchronization of WDC network. For 
this type of weighting scheme, the node weight intensity is 
s m = Kn E«eA m K = k nt a {K)m- According to our defi- 
nitions above, the coupling gradient between the pair of con- 
nected nodes, m and n, is, 
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For scale-free networks generated by the BA standard model 
(the model we have adopted for WDC network in Fig. Q3, the 
degrees of the nodes are not correlated, and we can roughly 
treat (k%) m as a constant, (fc«) m = (l/fe m ) £ fc« d- Us- 
ing this approximation, the above equation can be simplified 
into 



A#„ 



""m n 



h,a i,a k 1 + a i- 1 + Q 



k 



l+a 



km kn k m k^ 



(9) 



As shown by the above equation, the coupling gradient arises 
naturally when one of the conditions, a ^ —1 or k m ^ k n , 
is satisfied, which is generally met in weighted complex net- 
work. 

Separated by the value a = — 1, the function of the gra- 
dient at the two sides of the parameter space is completely 
different. Assuming that k m < k n in Eq. (0, in the region 
of a < —1, we statistically have s m > s n and, consequently, 
Ag mn < 0. Since s m oc k m , the gradient thus is pointing 
from the smaller-degree to larger-degree nodes. According to 
the result of Ref. II 1311 . this negative gradient will suppress 
synchronization and, moreover, in the negative regime with 
the increase of the gradient strength, (Ag mn ), the value of 
e c will be monotonically decreased. In a similar fashion, in 
the region of a > — 1, the gradient is pointing from the larger- 
degree to smaller-degree nodes, i.e. the gradient is positive. In 
this case, with the increase of the gradient strength, the value 
of e c will be also monotonically decreased. Combining the 
behaviors of e c in the two regimes, we thus predict a continu- 
ous decrease of the critical coupling strength as a function of 
the weight parameter, a, in the WDC network. 

The above predictions are verified by numerical simula- 
tions. By the scale-free network used in Fig. Q] we plot the 
variation of e c as a function of a in Fig. [3ja). Clearly, with 
the increase of a the network synchronization is found to be 
monotonically enhanced (a continuous decrease of e c ). In the 
inset of Fig. [3} a), we also plot the variation of the averaged 
network gradient, (Ag), as a function of a, where a contin- 
uous increase of (Ag) is evident. The enhanced synchro- 
nization by coupling gradient is more clearly shown in Fig. 
Efb), in which the value of e c is found to be monotonically 
decreased with the increase of the averaged network gradient. 

Our findings in WDC network, i.e., the results in Fig. [3] 
might give indications to the design of realistic networks 
where synchronization is importantly concerned, e.g., the 
world-wide airport network 113011 and the metabolic network 
lO . In a practical situation, due to the high cost of main- 
taining the weight gradient, a tradeoff between the network 
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FIG. 3: (Color online) For the same scale-free network as used in 
Fig. 1, under WDC weighting scheme, the variation of the critical 
coupling strength, e c , as a function of the weight parameter, a, in (a) 
and the average network gradient, (Ag), in (b). Inset: (Ag) versus 
a. Each data is averaged over 20 network runs. 



performance (synchronization) and the weight (coupling) gra- 
dient will be necessary. The results of Fig. [2 a) suggest that 
in the region of a < 2, the network synchronization can be 
significantly improved by a small gradient cost (small change 
of a); while in the region of a > 2, to achieve the same 
performance improvement, the gradient cost could be signif- 
icantly large. This might be one of the reasons why some 
realistic networks are found to be possessing weight parame- 
ter a « 0.5, e.g., the world-wide airport network 113011 and the 
E. coli metabolic network 1331. 



B. WDU-type Networks 

We next analyze the synchronization of WDU network. As 
indicated in Fig. |2j in this type of networks with the increase 
of the weight parameter, (3, the network synchronization may 
be either enhanced or suppressed. To have a close look to the 
relationship between e c on f3, we increase j3 from to 4 grad- 
ually, and checking the variation of e c . The numerical results 
carried on scale-free networks are plotted in Fig. Sl a )- In this 
figure, it is clearly shown that at about /3 C f=s 1.5, there ex- 
ists a maximum value for e c . It is worthwhile to note that the 
cap-shape variation of e c may have implication to the function 
of some realistic networks where network desynchronization 
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FIG. 4: (Color online) For the same scale-free network as used in 
Fig. 2(a), under WDU scheme, the variation of the critical coupling 
strength, e c , as a function of the weight parameter, /?, in (a) and the 
averaged network gradient, (|A<?|), in (b). Insets: (|A(?|) versus /3. 
Each data is averaged over 20 network runs. 



FIG. 5: (Color online) For the same small-world network as used in 
Fig. 2(b), under WDU scheme, the variation of the critical coupling 
strength, e c , as a function of the weight parameter, /?, in (a) and the 
average network gradient, (|A<?|), in (b). Insets: (|Ac?|) versus /3. 
Each data is averaged over 20 network runs. 



lT38h . instead of synchronization, is desired. Say, for instance, 
in the mobile communication network, simultaneous commu- 
nications may cause the information congestion and should be 
always avoided. An interesting empirical observation is that, 
in this weighted network, the weight of the network links fol- 
lows roughly a pow-law distribution with an exponent j3 « 1.9 
ll32tl . This value is quite close to (3 C in Fig. Ufa), which might 
indicate the possible application of weight-degree correlation 
in affecting network dynamics for this specific network. 

We now analyze the synchronization behavior of WDU net- 
work from the viewpoint of gradient network. It should be 
noted that, different to the case of WDC network, in WDU 
network the direction of the gradient is not uniquely defined. 
That is, the gradients on the network links may either point- 
ing from the larger-degree node to the smaller-degree node, 
or pointing in the opposite direction. As a matter of fact, 
in WDU network the gradient directions are determined by 
the weight details, instead of node degrees. That is, the posi- 
tive and negative gradients are coexisting in a WDU network. 
This arises a problem in using the averaged network gradient 
for the synchronization analysis, since for a large-scale WDU 
network, due to the random weight arrangement, we roughly 
have (Ag) ss 0. Regarding of this, in WDU network we use 

the averaged gradient amplitude, (|Ag|) = j^ 9 " 1 "^ , to 



characterize the gradient at the network level. With this adop- 
tion, our task now is shifted to exploring the relationship be- 
tween e c and (| Ag\). 

We first check the variation of (|Ag|) as a function of /3. 
The results are plotted in the inset of Fig. IS a). Different to 
the case of WDC network, it is observed that here the varia- 
tion of (| Ag\) has a cap-shape structure. (For WDC network, 
as a increases from 0, the value of (Ag) is monotonically in- 
creased, see Fig. [3fa).) In particular, at the critical parameter 
(3 C , we also find a maximum value for ( | Ag \ ) . The similar be- 
havior of e c and (| Ag\) seems to indicate the following: gra- 
dient deteriorates synchronization in WDU network. To check 
this out, we plot in Fig. |4jb) the dependence of e c on (|Ag|). 
Now it is clearly seen that, as the network gradient increases, 
the network synchronization is monotonically suppressed. 

The suppressed synchronization by gradient is also found 
in small-world networks under the WDU scheme. With the 
small-world network model used in Fig. [2] we calculate the 
variation of the critical coupling strength as a function of the 
weight parameter. The results are presented in Fig. 0a). 
Again, the maximum value of e c is found at some interme- 
diate value of (3. Similar to the case of scale-free network 
(Fig. |4j, the averaged gradient amplitude, (| Ag\), also shows 
a cap-shape structure (inset of Fig. |5ja)), with the maximum 
value of (| Ay |) appears at /3 C . The relationship between e c 
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and (|Ag|) is plotted in Fig. [5fb), where the suppression of 
synchronization by gradient is evident. 



C. Mechanism analysis 

Why synchronization is suppressed by gradient in WDU 
network, while is enhanced in WDC network? To answer this, 
we compare the structures of the gradient networks generated 
in the two types of networks, and analyzed their functions to 
the network synchronization. For a WDC network, according 
to the study of Ref. lf26ll . the gradient couplings of the net- 
work are organized into a spanned tree. More specifically, at 
the top of the gradient network there locates the largest-degree 
node of the network, which sending out gradients to its nearest 
neighbors (the second level). Each node at the second level is 
then sending out gradients to its own nearest neighbors (the 
third level), given that it has a larger degree to its neighbors. 
In this way, all the network nodes will be organized into a one- 
way coupled tree structure. An important feature of this gra- 
dient tree is that, starting from the rooting node (the largest- 
degree node), every node of the network can be reached by 
following the gradient couplings. Previous studies on gradi- 
ent network have shown that B25I1 . for such a spanned gradient 
tree, the network synchronizability will be monotonically in- 
creased with the increase of the gradient strength. From Eqs. 
dU and d6}, it is straightforward to see that in our model of 
WDC network, the gradient strength is increasing with the 
parameter a in a monotonic fashion. We thus attribute the 
enhanced synchronization in WDC network to the existence 
of a spanned gradient tree. 

To better describe the spanned gradient tree in WDC net- 
work, we adopt the simplified network model plotted in Fig. 
Oa), which might be also regarded as a motif of a large-scale 
network. According to our model (the Eq. ((3)), under WDC 
scheme the gradient direction will be pointing from larger- 
degree node to smaller-degree node on each link. Follow- 
ing this rule, we can construct the gradient network, which 
is schematically shown in Fig. |6|b). Please note that, to illus- 
trate, here we only show the gradient direction, while neglect- 
ing the gradient strength. In Fig. 0b), it is clearly seen that, 
leading by node 1, all the network nodes are connected by the 
gradient couplings in a hierarchical fashion. It is important 
to note that, while the exact gradient strength of each link is 
modified by the parameter a, the hierarchical structure of the 
spanned gradient tree is unchanged. In particular, in the gradi- 
ent network node 4 is always reachable by node 1, by passing 
the nodes 2 and 3. 

We go on to study the gradient couplings in WDU network. 
To manifest the feature of weight-degree non-correlation, and 
also for the sake of simplicity, we keep all other weights of the 
network links the same to those of Fig. |6jb) (WDC network), 
while changing only the weight between nodes 3 and 4 by a 
large increment, A. Taking into account this single change of 
the link weight, we now reanalyze the gradient couplings of 
the network. Since all other weights of the network are un- 
changed, all the gradients are the same to that of Fig. ISfb), 
except the one on the link between nodes 2 and 3. According 




(b) (c) 

FIG. 6: (Color online) A schematic plot on the distribution of the 
gradient couplings in WDC and WDU networks, (a) The simplified 
model of undirected, non-weighted network, (b) The distribution of 
the gradient couplings under the WDC-scheme. The network nodes 
and the gradient couplings are organized into a spanned gradient tree, 
in which each node is reachable from the rooting node (numbered 
1). (c) The distribution of the gradient couplings under the WDU- 
scheme, in which the weight on the link between nodes 3 and 4 is 
increased by A, while the remaining gradients keep identical to that 
of (b). As A exceeds some critical value, the gradient between nodes 

2 and 3 will switch its direction (the red arrow-line), manifesting a 
breaking of the gradient network. Dashed lines are links dominated 
by symmetric couplings (with negligible gradient). 

to our model of Eq. ([3J, the gradient from node 2 to node 3 is 
A(?23 = C23 (I/S3 — I/S2), with s the weight intensity defined 
in Sec. II. In WDC network (Fig. |6fb)), because S2 > S3, 
we thus have A523 > 0. This direction of gradient, however, 
will be switched in WDU network if A is large enough (Fig. 
13c)). In the WDU network, the gradient from node 2 to node 

3 is 5g' 23 = C23(l/sg — I/S2), with S3 = S3 + A. It is not 
difficult to find that, as long as A > S2 — S3, the gradient will 
be pointing from node 3 to node 2. That is, the gradient direc- 
tion is switched. Once switched, the gradient network will be 
separated into two sub-trees: one led by node 1 and the other 
one led by node 3. Previous studies have shown that l26ll . 
for such a breaking network, the network synchronizability 
will be monotonically deteriorated by increasing the gradient 
strength. Based on this analysis, we attribute the deteriorated 
synchronization in WDU networks to the breaking of the gra- 
dient network. 



V. DISCUSSION AND CONCLUSION 

The main purpose of the present paper is to highlight the 
important role of the weight-degree correlation played on the 
network synchronization. While synchronization in weighted 
complex network has been extensively studied in literature, 
to the best of our knowledge, it is the first time that the cor- 
relation between the link weight and node degree be investi- 
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gated. To demonstrate the effect of weight-degree correlation 
on the network dynamics, we have employed the generalized 
Kuramoto model as the platform and, following the tradition, 
adopting the scheme of normalized coupling strength. The 
main findings of the paper, however, are general and are ex- 
pected to be observable for other network models and cou- 
pling schemes as well, say, for instance, the complete network 
synchronization. 

Motivated by the empirical observations, we have adopted 
two different weighting schemes for the WDC and WDU net- 
works. The weight parameters, a and f3, thus have different 
meanings and, in general, should not be directly compared. 
For this regard, the comparison we have made in the present 
work is actually on the different dynamical responses of the 
two types of networks to the change of the weight hetero- 
geneity, which should not be understood as a comparison of 
the network synchronizability. Direct comparison of the two 
weighting schemes might be possible for the BA-type scale- 
free networks, in which we roughly have j3 = —7/ (2a). 
However, to keep this relation, in numerical simulations the 
network size should be very large (larger than 1 x 10 6 ), which 
is out of our current computational ability. 

While gradient coupling has been identified as impor- 
tant for network synchronization, previous studies have been 
mainly concentrated on the type of WDC network ll28ll . As 
such, a general observation is that, by increasing the gradient 
strength, the network synchronizability is monotonically in- 
creased. This observation, as we show in the present paper, is 
invalid for the case of WDU network. In particular, in WDU 
network we find that, instead of enhancement, the increase 
of the gradient could actually suppress synchronization under 
some conditions. This finding could be helpful to the design 



of realistic networks which, depending on the network func- 
tionality, may either in favor of network synchronization or 
network desynchronization. Now, besides the aspect of link 
weight, to manipulate the network synchronization, we can 
also adjust the correlation between the link weight and the 
node degree. 

We wish to note that, despite of the recent progresses in the 
onset of synchronization of networked phase oscillators, at the 
present stage we are still short of a theory to predict precisely 
the point of onset of synchronization ifXOll - This is specially 
the case when the link weight is of heterogeneous distribution, 
e.g. the WDC and WDU networks discussed here. While our 
main findings are based on numerical simulations, it is our 
believing that these findings are general and stand for other 
network models and dynamics as well. 

In summary, we have investigated the onset of synchro- 
nization in weighted complex networks, and explored the in- 
fluence of weight-degree correlation on the network synchro- 
nization. An interesting finding is that, as the weight parame- 
ter increases, in WDC network the synchronization is mono- 
tonically enhanced, while in WDU network the synchroniza- 
tion could be either enhanced or deteriorated. This finding 
is obtained from extensive simulations, and is explained from 
the viewpoint of gradient network. Our study shows that, be- 
sides the link weight, the weight-degree correlation is also an 
important concern of the network dynamics. 
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